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ABSTRACT 

Fisher Information Matrix methods are commonly used in cosmology to estimate the 
accuracy that cosmological parameters can be measured with a given experiment, 
and to optimise the design of experiments. However, the standard approach usually 
assumes both data and parameter estimates are Gaussian-distributed. Further, for 
survey forecasts and optimisation it is usually assumed the power-spectra covariancc 
matrix is diagonal in Fourier-space. But in the low-redshift Universe, non-linear mode- 
coupling will tend to correlate small-scale power, moving information from lower to 
higher-order moments of the held. This movement of information will change the 
predictions of cosmological parameter accuracy. In this paper we quantify this loss of 
information by comparing naive Gaussian Fisher matrix forecasts with a Maximum 
Likelihood parameter estimation analysis of a suite of mock weak lensing catalogues 
derived from N-body simulations, based on the SUNGLASS pipeline, for a 2-D and 
tomographic shear analysis of a Euclid-like survey. In both cases we find that the 68% 
confidence area of the VL m — a% plane increases by a factor 5. However, the marginal 
errors increase by just 20 to 40%. We propose a new method to model the effects 
of nonlinear shear-power mode-coupling in the Fisher Matrix by approximating the 
shear-power distribution as a multivariate Gaussian with a covariance matrix derived 
from the mock weak lensing survey. We find that this approximation can reproduce 
the 68% confidence regions of the full Maximum Likelihood analysis in the il m — o~8 
plane to high accuracy for both 2-D and tomographic weak lensing surveys. Finally, 
we perform a multi-parameter analysis of f2 m , erg, h, n s , wq and w a to compare the 
Gaussian and non-linear mode-coupled Fisher matrix contours. The 6-D volume of 
the la error contours for the non-linear Fisher analysis is a factor of 3 larger than 
for the Gaussian case, and the shape of the 68% confidence volume is modified. We 
propose that future Fisher Matrix estimates of cosmological parameter accuracies 
should include mode-coupling effects. 
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1 INTRODUCTION 

Measuring cosmological parameters and the equation of 
state of dark energy to high accuracy is thegoal of many cur- 
rent and upcoming experiments (e.g. PSIq DEjQ, Euclid^, 
HALC0, LSSuf] and WFIRSlQ). Variables such as the size 
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and depth of the survey (amongst other things) have a sig- 
nificant effect on the ability of a survey to constrain cosmo- 
logical parameters. Consequently, significant effort must be 
spent in accurately predicting what these telescopes will see, 
well before construction begins. This will allow us to both 
to influence the design phase and to understand the capa- 
bilities of the instrument once the design has been set. In 
order to make predictions for these upcoming missions, sta- 
tistical tools must be used to estimate the accuracy they will 
be able to achieve. The current standard for prediction uses 
the F isher matrix methodology (|Tegmark et al.l ll997: Fisher 
Il935t) . Traditionally, Fisher matrices have been generated 
with data covariance matrices that assume an underlying 
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Gaus s ian matter a nd radiation distribution jTegmark et al.l 
1 19971 ; lKnoxlll997h , which is accurate for CMB estimates 
when the Universe was still linear. However, this is not an ac- 
curate representation of the low-redshift Universe at smaller 
scales. Error bars generated using these Gaussian assump- 
tions may be biased when compared with those generated 
using methods that account for the non-Gaussian, non-linear 
nature of the Universe. 

In this paper we compare the la, two-parameter £l m —as 
error estimates from the full maximum likelihood analy- 
ses with a Gaussian Fisher analysis. Using weak gravi- 
tational lensing power spectrum analysis of mock galaxy 
shear catalogues gen erated with the SUNGLASS pipeline 
(jKiessling et al.ll201ll ). we show the importance of using ac- 
curate non-linear covariance matrices when estimating er- 
rors for future experiments. Analytic approximations of the 
correlation function covariance matrix under these assump- 
tions have been shown to underestimate the errors on cosmic 
variance by a factor of up to ~ 30, w hich makes breaking th e 
flm — as degeneracy more difficult (|Semboloni et al.l l2007). 
Simulations are able to provide accurate covariance matrices 
because they do not make assumptions about the underlying 
Gaussianity of the Universe and consequently, the resulting 
covariance matrices include the off-diagonal components. 

Fisher matrix analyses are attractive because of their 
relative speed and minimal computational requirements 
when compared with maximum likelihood estimates. How- 
ever, this is offset by the loss of information due to Gaus- 
sian assumptions in the generation of the Fisher matrix. 
To compensate for some of these Gaussian assumptions, we 
propose using a non-Gaussian weak lensing shear covariance 
matrix generated from simulations to generate the Fisher 
matrix. The resulting Fisher matrix still assumes a multi- 
variate Gaussian parameter estimate distribution and power 
spectrum distribution. However, it now contains the non- 
linear information found in the off-diagonal components of 
the covariance matrix from the simulations, giving us a 'non- 
linear' Fisher matrix. In this paper, we will show the effect 
of this simple modification to the calculation of the Fisher 
matrix on the error estimates and compare this with the 
maximum likelihood and Gaussian Fisher error contours in 
the Q. m — as plane. We also perform a multi-parameter anal- 
ysis of fl m , as, h, n s , »o and w a and compare the Gaussian 
and non-linear Fisher matrix error estimates. 

It is possible to produce data vectors for covariance 
matrices by performing a simple 2-D binning of the galax- 
ies in the survey. However, further information may be 
gained by splitting the distribution up in to redshift bins 
and performing a to mographic analysis (e.g. Hu 1999, 2002; 
I Jain fc Taviorll2003r i. In this work, we perform both 2-D and 
3-bin tomographic analyses in order to generate covariance 
matrices that are used to calculate the maximum likelihood 
estimates and the Gaussian and non-linear Fisher matrices. 

The outline of this paper is as follows. Section [2] will de- 
tail how the simulations and mock galaxy catalogues were 
generated. Section[3]will introduce the Fisher matrix formal- 
ism and the maximum likelihood formalism. Section [4] shows 
the results of the analyses on the maximum likelihood esti- 
mates and the Gaussian and non-linear Fisher matrices from 
the 2-D analyses in Section T4. II and the tomographic analy- 
ses in Section r4.2l Section ^. 41 compares the multi-parameter 



N Area (sq deg) z max n g /sq arcmin z med 



100 



100 



1.5 



15 



0.82 



Table 1. Table of parameters for the mock galaxy shear cata- 
logues used in this paper. N is the number of independent light- 
cones, z max is the maximum redshift in the lightcone and n g /sq 
arcmin is the number of 'galaxies' per square arcmin in the cata- 
logue and z me< j is the median redshift of the catalogue. The suite 
of lightconcs is used together to form a survey with an effective 
area of 10,000 sq. deg. 



Gaussian and non-linear Fisher estimates. Finally, a sum- 
mary of the findings will be presented in Section [S] 



2 DETAILS OF THE SIMULATIONS 

The suite of weak lensing simulations used in this work was 
generated usin g the SUNGLASS pip eline (for a detailed in- 
troduction, see iKiessling et~al]|201lh . We have 100 indepen- 
dent simulations generated w ith the cosmol ogical N-body 
simulations code GADGET2 l|Springell [20051 ). The simula- 
tions were made with a flat concordance A CDM cosmol- 
ogy, c onsistent with the WMAP 7- year results l|jarosik et al.l 
1201 lb : Q m = 0.272, Q A = 0.728, fl b = 0.045, as = 0.809, 
n s = 0.963 and h = 0.71 in units of 100 km s" 1 Mpc" 1 . 
There are 512 3 particles in a box of 512/i _1 Mpc which leads 
to a particle mass of 7.5 x 1O 1O M . The simulations were all 
started from a redshift of z — 60 and allowed to evolve 
to the present with 26 snapshots being stored in redshifts 
0.0 < z < 1.5. 

Lightcones were generated through the simulation snap- 
shots to determine the average convergence in an angular 
pixel using the 'no r adial binning' method introduced in 
IKiessling etaTl (|201ll ): 



K(r k ,r s 



Aft 

k 



P n(r k )r% 



dr K(r, r 3 



(1) 



where n is the number density of particles in the simulation, 
AQ P = A9 x A0 y is the pixel area, r s is the comoving radial 
distance of the lensing source plane and rk is the comoving 
radial distance of each particle, k, in the lightcone. K(r,r s ) 
is the scaled lensing kernel: 

3HoQ m (r s - r)r 



K(r, r s 



2c 2 



(2) 



r a a{r) 

The convergences were calculated using 2048 2 azimuthal 
bins on source redshift planes that were separated by z = 0.1 
to create 15 planes from 0.0 < z < 1.5. Once the con- 
vergences were calculated, shear values were determined in 
Fourier space, where the shear is given by 7 = 71 + ij2- 

The shear and convergence values in these source red- 
shift planes were interpolated back onto the individual par- 
ticles in the lightcones to generate mock galaxy shear cata- 
logues. The B-modes (the unphysical, imaginary component 
of the convergence) in the mock catalogues were calculated 
directly from the shear. These catalogues were then sam- 
pled to_j^p_roohic^_a staxidard galaxy redshift distribution 
(e.g. lEfstathiou et al.lll99ll 'l. 
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n(z) oc z a exp 



(-) 



(3) 



where z is the redshift, Zo, a and 77 set the depth, low-redshift 
slope and high-redshift cut-off for a given galaxy survey. We 
take a = r\ = 2 and 20 = 0.78 in this work, which gives 
a median redshift z me d = 0.82. It is assumed that galax- 
ies trace the dark matter distribution perfectly and the fi- 
nal mock galaxy shear catalogues contain 15 galaxies per 
square arcminute. There is no ellipticity noise in the cata- 
logues, however there is a shot-noise contribution related to 
the discrete sampling of the particles in the mock catalogues. 
Table [1] summarises the mock galaxy catalogues used in this 
work. 



3 METHODOLOGY 

The shear, 7, convergence, k, and B-mode, /3, fields are re- 
lated to each other in Fourier-space on a fiat-sky by 



(4) 



where ipi is the angle between the angular wave-vector and 
an axis on the sky. For each of the mock galaxy shear cata- 
logues the shear, convergence, and B-mode auto- and cross- 
power spectra have been estimated. The tomographic weak 
lensing shear, convergence and B-mode power cross-spectra, 
for two different source redshifts z and z' , are given by 

C'P{z,z') = ( 7l (^) 7l (^')> + < 72 (£,*) 72 (5) 

cr (*,**) = w,z)K{t,z')), (6) 

Cf(z,z') = {I3(£,z)/3(£,z')). (7) 

The power spectra are related to each other by C~p (2, z') = 
C( K (z,z') + Ci (z, z'). The auto-spectra are calculated 
when z' = 2. In practice we bin the sources into redshift 
slices. For a survey that has N z slices the expectation- value 
of the tomographic shear power spectrum in redshift bins 
labeled i and j, is given by 



dr 



^,r)m{r) gj {r), (8) 



3 4c 4 J Q a 2 (r) 

i|Kaiserl Il992l ; | Joachimi fe Schneider! [2009). where r is the 
comoving distance, r max is the maximum comoving distance 
and P(£/r, r) is the 3D matter density power spectrum. The 
lensing efficiency function is 



9i(r) = / dr Pi (r') (l - ~) 



(9) 



where pi is the normalised probability distribution of 'galax- 
ies' in the bin. 

When determining the angle averaged shear power spec- 
trum in the simulations, we must take into consideration the 
conventions used in the Fourier transform software FFT"wQ. 
Thus, the discretised tomographic shear power spectrum be- 
comes 



l(l + l)Cff(£) 
2tt 



E 71 (I, gi)7l (£, Z J ) + 72 (I, Zi) 7 2 (I, Zj ) 
n 2 Aln^ ,{ > 



7 The Fastest Fourier Transform in the West 
http://www.fftw.org 



where CJ 1 is the estimated power, rib is the total number of 
bins in the Fourier transform, 2; is the i th redshift slice, and 
Aln^ is the thickness of a shell in log ^-space. The modes 
in this power spectrum are corrected for mode discreteness 
errors by scaling by the expected number of modes. To com- 
pactify the notation we shall denote the tomographic shear 
power by C]^ (£), where X = is a pair of redshift slices. 
For our 2-D analysis there is only one bin so i = j. Using 
the information from these power spectra, we have the the 
information required to make cosmological parameter esti- 
mates. 



3.1 Shear power covariance matrix generation 



The shear cross-spectra covariance matrix, M^, x is defined 
by 



(11) 



where ACJ7(£) = CJP '(£) - {C]? '(£)} and angled-brackets 
denotes ensemble averaging. The covariance matrix is an im- 
portant element of parameter estimation, containing infor- 
mation on the strength of the correlations between variates, 
in this case the shear power spectrum modes. The accuracy 
of this matrix improves by increasing the numbers of real- 
isations included in the calculation of the matrix. For this 
work, 100 independent realisations of 2-D power spectra pro- 
vides an accuracy of AC'e ~ 10% (as shown in Figure 7 of 
iKiessling et"alll201lh . 



3.1.1 Gaussian covariance matrix 

If the shear field is isotropic and we assume the shear field 
is Gaussian we can write down an equation for the tomo- 
graphic shear power covariance matrix and make a correc- 
tion for the fraction of the sky covered by the survey: 



(2£ + l)/ sky 



([C^ M + N^S^lCJ^i) + Nj (£)6 



+ [(%» {£) + Ni (£)5$ ] [CTS (I) + Nj [£)5f k ] , (12) 



Ml]* 1 



where / s k y is the fraction of the sky covered by the survey 
and Ni is a shot-noise term due to intrinsic ellipticity in the 
shear field. For the purposes of this work although we have 
discrete galaxy density, we set the intrinsic shear variance to 
zero so that Ne ~ 0. In principle our results can be combined 
with a Gaussian noise covariance to model different weak 
lensing surveys. 



3.1.2 Simulation covariance matrix estimation 

The covariance matrix of the tomographic shear power spec- 
trum for the suite of mock galaxy shear catalogues is esti- 
mated by 



Ale, 



N-l^> 

N 



ACT?{£)AC]?,{£') 



(13) 



where N is the number of mock catalogue realisations and 
here AC]7(l) = d%> '(£) - where {CJ7 {£)} is the 

ensemble average of the tomographic shear power spectrum 
across all realisations. 
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We calculate inverse covariance matrix using a singu- 
lar value decomp osition (SVD) on the covariance matrix 
l|Press et al.lll992l ). However, the resulting inverse is biased 
due to n oise. To correct for t his, we multiply the inverse by 
a factor l|Hartlap et al.ll2007l ): 



N s -N p 



Ns _ l i , (14) 

where Ns is the number of realisations, N p is the total num- 
ber of bins in all of the power spectra and [M x , x ] _1 is an 
unbiased estimation of the inverse covariance matrix. 

3.2 Fisher matrix generation 



The Fisher matrix is given by 



Fi 



EE 



dc%" 



dcrg 



(15) 



where i and j label cosmological parameters (e.g. Q m and 
as) and the partial derivatives are the gradient of the expec- 
tation value of the shear power spectrum in parameter space. 
In order to determine tha t gradients of the power spect ra, we 
use a five-point function l|Abramowitz fc Steeun|[l968l . equa- 
tion 25.3.6) and the ensemble average power spectra. The 
theoretical prediction for the power spectrum was gener- 
at ed using a code provided by Benjamin Joachimi (as tested 
in iKiessling et all 1201 ll; Ijoachimi fc Schneiden 120091 ) . This 
code uses the lSmith et al.l (120031') non -linear power spectrum 
estimate, the Eisenstein fc Hul(ll998l ) matter transfer func- 
tion and a numerically calculated linear growth factor. All 
further theoretical power spectrum predictions in this paper 
were performed using this same code. 

The Fisher matrix provides parameter estimates 
through two methods. The first by calculating the area of 
the Fisher matrix error ellipse that encloses a two-parameter 
68% confidence limit in the two parameter plane. The in- 
verse of this area is proportiona l to the Figure-of-Merit that 
is often quoted in studies (e.g. A 



brecht et al.ll2006l: IWanel 



12008! : lAlbrecht et al~ll2009l ; ISlosarll2010l ~ The second is the 
single-parameter marginal errors, which are given by 

1 1/2 



AO, 



(16) 



3.3 Maximum likelihood parameter estimation 

To perform a maximum likelihood analysis on the suite of 
mock shear catalogues, we use a Gaussian likelihood estima- 
tor. Despite the fact that the simulations are non-Gaussian, 
using a Gaussi an likelihood has bee n shown to produce ac- 
curate results (|Kiessling et al.l[201ll ) . The likelihood is given 
by 



L(Cr\a 8 ,n m ) 



(2n) N / 2 (detM xx ') 1 / 2 



exp 



' 2 



(17) 



where N is the number of independent mock catalogue re- 
alisations, and 



M«' X ' 



AC^(£'), 



(18) 



where AC* 77 (£) = C 77 '(e) - (C 77 00), and (C, 77 ) is the ex- 
pected angular power spectrum given by equation (8). While 




100 1000 

Wovenumber I 

Figure 1. Mean 2-D shear power spectrum for the suite of un- 
binned mock galaxy shear catalogues. The smooth (red) line is 
the theoretical prediction for the power spectrum, the (black) 
line shows the mean power spectrum for the suite of mock cata- 
logues with errors on the mean shown. The diagonal (dark blue) 
dashed line shows the shot noise prediction and the (magenta) 
triangles show the measured B-modes. The bottom panel shows 
the the percentage difference of the data from the expected power 
spectrum with errors. 



this likelihood analysis is computationally expensive, it takes 
into account the full non- Gaussian and non- linear nature of 
the simulations and should provide the most accurate error 
estimates. The maximum likelihood is related to the Fisher 
matrix by 



f 



d 2 InL 



(19) 



4 RESULTS 

4.1 2-D two-parameter analysis 

The first step to making parameter estimates from the mock 
galaxy shear catalogues is to determine the 2-D shear angu- 
lar power spectra for each realisation. Figure [1] shows the 
mean power spectra from the suite with the (red) long line 
showing the theoretical prediction, the (black) line showing 
the mean measured power spectrum from the suite of cat- 
alogues with errors on the mean. The (magenta) triangles 
show the measured B-modes and the dashed (dark blue) line 
shows the shot-noise arising from discrete particle sampling 
in the mock catalogues. The bottom panel shows the per- 
centage difference of the data from the expected power spec- 
trum. The simulations recover the expected power spectrum 
within 5% between 150 < I < 2000. 

The shot noise prediction for these catalogues was de- 
termined by filling a suite of simulation volumes with ran- 
domly placed particles to mimic Poisson noise. The SUN- 
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500 1 ooo 

Wavenumber i 



2000 



Figure 2. Diagonal components of the covariance matrix for the 
simulations (black line) and a Gaussian field with the same power 
spectrum (red dashed line). 



GLASS pipeline was run on these noise simulations and the 
power spectrum of the mock catalogues is the dashed line in 
the figure. We know that the shot-noise in these catalogues 
is not purely Poissonian becaus e the simulat ions start from 
a glass pre-initial distribution |Whit3[l994). However, the 
noise tends toward Poissonian as the simulations evolve. We 
are not able to model this noise ac curately because it evolves 
with structures in the simulation dBaugh et al.|[l995l ). Thus, 
this Poisson estimate should be a reasonable approximation 
for the noise in these catalogues. 

A data-vector was constructed from the shear band- 
power spectra with wavenumbers from 150 < I < 1500. 
From this, we estimated the covariance matrix in equation 
(|13p . Similarly, a Gaussian covariance matrix was generated 
using equation (fl2)) . 

To test the effect of the non-Gaussian nature of the 
simulations covariance matrix, we looked at the diagonal 
components of both the simulation and Gaussian field co- 
variance matrices and determined their ACe values. For the 
case of the simulations this is simply the square root of the 
diagonal components. For the Gaussian field this is 



ACe 



V2Ct 



^(2£+l)/ sky Aln<? 



(20) 



Figure shows the diagonals of the covariance matrices 
as a function of wavenumber I. The (red) dashed line is from 
the Gaussian error and the (black) line is from the simula- 
tions. At low wavenumbers, which is still in the reasonably 
linear regime, the two errors agree reasonably well. However, 
in the higher wavenumber, non-linear regime, there is a fac- 
tor of three difference between the errors, suggesting that 
the non- Gaussian contribution to the covariance is signifi- 
cant. 

Fisher matrices in the Sl m — as plane were calculated 




0.26 



0.27 



0.28 



0.29 



Figure 3. Comparison of the la two-parameter contours for the 
— erg parameters for Gaussian Fisher contour (inner orange 
shaded ellipse) with the non-linear Fisher contour (middle grey 
shaded ellipse) and the full maximum likelihood analysis (thick 
black line contour). The (red) point in the middle of the ellipses 
represents the fiducial Q m — as parameters in this calculation. 



using equation (|15[1 for both the Gaussian and simulation co- 
variance matrices. The Fisher matrices were multiplied by 
100, the total number of independent mock catalogues, to 
provide an error estimate for a survey of 10,000 square de- 
grees. Additionally, the maximum likelihood was calculated 
using the simulation covariance matrix and the likelihoods 
from the 100 mock c atalogues were combined (as shown in 
iKiessling et al.ll201ll ). For all of these calculations it is as- 
sumed that all other cosmological parameters are known. 

In Figure [3] the two-parameter la contours for the 
fim — as plane are shown. The inner (orange shaded) el- 
lipse shows the Gaussian Fisher estimate, the outer (grey 
shaded) ellipse shows the non-linear Fisher estimate and the 
thick (black) line contour is the combined maximum likeli- 
hood estimate. The red point in the middle represents the 
fiducial fl m and as parameters. The size and shape of the 
simulation Fisher ellipse and the maximum likelihood con- 
tour are very similar, showing that the non-linear Fisher 
calculation is a good method for making future cosmologi- 
cal parameter estimates, provided the off-diagonal terms are 
included. The area of the non- linear Fisher estimate is 5.1 
times larger than the area of the Gaussian Fisher estimate 
which is a clear indication that the off-diagonal terms in the 
covariance matrix are essential for accurate estimation of 
the parameter errors. 

In addition to the area of the contours we also deter- 
mined the marginal errors on the parameters from the Fisher 
matrices. The non-linear tt m marginal error is 1.23 times 
larger than the Gaussian marginal errors and as is 1.17 times 
larger. From Figure [2] we would expect the marginal errors 
to be much larger. However, the off-diagonal terms in the 
simulation covariance matrix act to reduce the marginal er- 
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Figure 4. Marginal errors of Q m and erg as a function of maxi- 
mum wavcnumber in the Fisher calculation. The dashed lines are 
the marginal errors on f2 m as a function of maximum wavenumber 
J maI and the continuous lines are for erg. 

rors but increase the total area. This is demonstrated by 
setting the off-diagonal terms in the simulation covariance 
matrix are to zero. In this case the non-linear Fisher ellipse 
becomes far narrower but the marginal errors increase sig- 
nificantly. 

Figure [3] shows how the marginal errors change in 
the non-linear Fisher matrix as a function of maximum 
wavenumber f max . The marginal errors for both the non- 
linear and Gaussian Fisher matrices are remarkably similar 
across all values of i" m ax with the largest gain of information 
occurring for both f2 m and erg between 250 < ^ ma x < 500 
and again at ^ max > 1000. This similarity is due to an effect 
that the off-diagonal components of the simulation covari- 
ance matrix are having on the error ellipse, as discussed 
earlier. This figure shows that if the marginal error is the 
value of interest, the Gaussian Fisher matrix appears to 
provide a result that is comparable with the error obtained 
with the non-linear Fisher error estimate at wavenumbers 
between 150 < ^ max < 1500. However, when the errors are 
marginalised over two parameters, the area of the contours 
shows that the errors are being underestimated by a factor 
of five. 



4.2 Tomographic shear power 

It is possible to perform the same analysis that was per- 
formed in Section T4. II on tomographically-binned mock cat- 
alogues. Tomography introduces extra information into the 
analysis. In this analysis the suite of mock galaxy shear cata- 
logues was split into three redshift bins; bin 1: 0.0 ^ z < 0.5, 
bin 2: 0.5 < z 1.0 and bin 3: 1.0 < z < 1.5. 

Figure[5]shows the auto- and cross-power spectra for the 
3-bin tomographic analysis of the mock galaxy shear cata- 



logues. If we first focus on the auto-power spectra (Cn, C22 
and C33), we can see that the wavenumbers that we model 
accurately increase with redshift. For Cn, we find agree- 
ment for wavenumbers from 150 < £ < 700. The power 
spectra from C22 and C33 are slightly higher than the ex- 
pected power spectrum but still within 3% and are reliable 
up to £ = 1000 before shot-noise becomes dominant. 

The cross-bins C12, C13 and C23 are damped which is 
due to the shot noise dominance in the lower redshift bin. 
Consequently, we only use the results of these power spectra 
up to the I range accurately recovered by the auto-power 
spectra of the lower redshift bin. 

The power spectra in the auto- and cross-bins were 
turned into a data-vector with wavenumbers from 150 < 
I < 1800 included. A data covariance matrix was generated 
and from this, the correlation coefficient matrix where the 
correlation coefficients are given by 

rff = . BEL . (21) 

This matrix shows how (anti-)correlated each of the £ modes 
are. 

Figure [6] shows the correlation coefficient matrix for 
the tomographic data vector. The 'block' nature to the 
matrix indicates each tomographic bin pair. Each of the 
blocks along the diagonal represents the auto-correlation 
between each of the tomographic power spectrum analy- 
ses. The off-diagonal blocks show the cross-correlations be- 
tween the tomographic pair power spectra. As expected, the 
higher wavenumbers in each block are highly correlated. Ad- 
ditionally, the Cn auto-power spectrum is highly correlated. 
The C22 auto-power spectrum is significantly less correlated 
and the C33 auto-power spectrum bin has very low corre- 
lations between the £ modes. The indication of this is that 
the lower redshift tomographic bins should contribute fewer 
power spectrum bins to a covariance matrix being used for 
any kind of analysis. 

With the information gained from these power spectra 
and the correlation coefficient matrix, we are able to assem- 
ble a data-vector that contains the Ce from wavenumbers 
that are accurately reproducing the expected power spec- 
trum and have a reasonably low correlation. Thus, we select 
30 bins in total from the three auto-power spectra, with 
the highest wavenumber being around £ = 1000 and just 10 
bins in total from the cross-power spectra, with the high- 
est wavenumber being around £ = 500. A covariance ma- 
trix was generated using this data- vector. A full maximum 
likelihood analysis was performed with this covariance ma- 
trix and a non-linear Fisher matrix was also calculated. An 
equivalent Gaussian covariance matrix was also generated 
and the Gaussian Fisher matrix determined. 



4.3 Two-parameter tomographic analysis 

Figure [7] shows the lcr, two-parameter contours for the Gaus- 
sian Fisher (orange shaded ellipse) and the non-linear Fisher 
(grey shaded ellipse). The thick black line shows the simu- 
lation maximum likelihood analysis. As in the case for the 
2-D analysis, the non-linear Fisher contour is very close to 
the maximum likelihood contour. The area of the non-linear 
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Figure 5. 2-D shear power spectra for a 3-bin tomography analysis of the mock galaxy shear catalogues. The long (red) line is the 
theoretical prediction for the shear power spectrum, the (black) line is the mean power spectrum for the 100 mock catalogues with errors 
on the mean, the (dark blue) dashed line is the shot noise estimate and the (magenta) triangles are the measured B-modes. 
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Figure 6. Correlation coefficient matrix for 3 tomographic bins. The matrix is arranged so that each sub-square is a tomographic bin 
pair. 



Fisher contour is 5.1 times larger than the area of the Gaus- 
sian Fisher contour. The marginal error of the non-linear 
Fisher estimate is 1.39 times larger in £l m and 1.24 times 
larger in <rg. These numbers are very similar to those calcu- 
lated in the 2-D analysis. 

The area of the contours and the size of the marginal 
errors in the tomographic Fisher and likelihood estimates is 
smaller than those in the unbinned analyses, showing that 
the tomographic analysis does indeed contain more informa- 
tion. 



4.4 Multi-parameter Fisher analysis 

In the previous section, we showed that a Fisher matrix 
analysis in the fl m — as plane accurately determines the er- 



rors on the parameters when a non-Gaussian data covariance 
matrix is used in the calculation. This was shown by com- 
paring the 68% error contour with a maximum likelihood 
analysis of the simulation suite. In this Section we perform 
a Fisher matrix analysis over the cosmological parameters: 
f2 m , erg, h, n s , wo and w a - Note that Q,a = 1— fi m - We com- 
pare the Gaussian Fisher analysis with the non-Gaussian 
(simulation) Fisher analysis. Based on the findings in the 
two-parameter analysis, we assume that the Fisher errors 
calculated with the simulation data covariance matrix are 
consistent with a maximum likelihood analysis of the simu- 
lations over the same multi-parameter space. 

For a Euclid-like survey (Table [1]), it makes little sense 
to perform a 2-D analysis over multiple parameters. We 
found that the size of the marginal errors was so large as 
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Figure 7. Comparison of the lcr, two-parameter contours for the 
— o"8 parameters for the Gaussian Fisher contour (inner or- 
ange shaded ellipse), the simulation Fisher contour (middle grey 
shaded ellipse) and the full simulation maximum likelihood anal- 
ysis (thick black line contour). The (red) point in the middle of 
the contours represents the fiducial Q m — parameters in this 
calculation. 
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0.107 
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Table 2. Marginal errors of a 3-bin tomographic, six-paramctcr 
Fisher matrix analysis using both Gaussian and non-linear Fisher 
matrices. For this configuration, the Gaussian marginal errors are 
always larger than the non-linear Fisher errors. However, the non- 
linear Fisher volume is 3.7 times larger than the Gaussian volume. 

to provide no constraining information (e.g. Awq = 11.3 
and Aw a = 49.8). However, the additional information pro- 
vided by performing a tomographic analysis of the power 
spectra yielded far better constraints and illuminated some 
features of the Fisher matrix analysis that were not immedi- 
ately obvious. Figure [8] shows the projected two-parameter, 
lcr contours, marginalised over the the multiple parameters 
from the Fisher matrix analysis. The (red) point at the cen- 
ter shows the fiducial parameters used in the analysis, the 
(grey) shaded ellipses show the non-linear Fisher analysis 
and the dashed (orange) ellipses show the Gaussian Fisher 
analysis. In all cases, except the Q, m — as plane, the projected 
area of the Gaussian contours is larger than the area of the 
non-linear Fisher contours which gives the impression that 



the error estimates from the Gaussian Fisher analysis are 
more conservative the non-linear estimates. The marginal 
errors of the Gaussian contours are larger in all cases (Table 
[2J. However, volume of the six-parameter space generated, 

Y oc v / det(F- 1 ), (22) 

shows that the non-linear Fisher volume is 3.7 times larger 
than the Gaussian volume, which is what we would expect 
given the off-diagonal terms included in the simulation data 
covariance matrix. The implication of this is that although 
the projected areas of the Gaussian Fisher appear larger, 
the overall volume is smaller. This can be explained with an 
example (in three dimensions for simplicity of explanation): 
Take a spherical ball and a thin plate with a slightly larger 
radius than the ball. If we look at the ball from any of the 
three axes, the area will appear to be the same circle in each 
projection. If we place the plate at a 45° angle to each of the 
three axes and then look at it in projection, it will appear to 
have a larger area than the ball in all projections, however 
we know that the volume of the plate is far smaller than 
the ball. With this knowledge in hand, it is easy to see how 
misleading the projections of the Fisher matrices can be in 
these complex multiple-parameter spaces. 

The eigenvalues and eigenvectors of the Fisher matrices 
were determined to find the equation for the plane with the 
most information in the full six-parameter space. For the 
Gaussian, this is: 

X = 0.865 Q. m + 0.498 a s + 0.049 n s + 0.030 h - 

0.035 wo - 0.006 w a , (23) 

and for the simulations: 

Y = 0.841 fi m + 0.536 a 8 + 0.034 n s + 0.028 h - 

0.056 w -0.015 w a . (24) 
When using the fiducial parameters, 

X = 0.741, (25) 

Y = 0.771, (26) 

and the errors on the planes are AX — 8.6 x 10 -8 and AY" = 
1.4 x 10 -6 . This shows an order-of- magnitude difference in 
the thickness of the 6-D error ellipses. 

5 DISCUSSION AND CONCLUSIONS 

This paper makes a comparison between Fisher matrices 
in the O m — as plane and six-parameter Fisher matrices 
(f2 m , as, h, n s , wo and w a ), generated using covariance 
matrices from a Gaussian random field and from full N- 
body simulations. The la two-parameter contours from the 
two-parameter Q. m — as Fisher matrices are also compared 
with the contour from a maximum likelihood analysis of the 
full suite of simulations. 

This work uses the SUNGLASS pipeline to generate 100 
independent simulations of 512/i -1 Mpc with 512 3 particles 
and a standard ACDM cosmology. The pipeline turns these 
simulations into 100 independent mock galaxy shear cata- 
logues of 100 square degrees and a galaxy redshift distribu- 
tion with a median of z m = 0.82 and 15 'galaxies' per square 
arcminute. When these catalogues are combined, they pro- 
vide an effective survey area of 10,000 sq. deg. 
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Figure 8. Multiple parameter Fisher analysis. The (grey) shaded ellipse is the non-linear Fisher contour and the (orange) dashed ellipse 
is the Fisher contour for the Gaussian Fisher matrix. The volume of the 6-dimensional space is larger in the simulation Fisher matrix, 
however some of the projected contours appear larger for the Gaussian Fisher matrix. This could be very misleading. 



We perform both 2-D and 3-bin tomographic angular 
shear power spectrum analyses on each of the mock cata- 
logues and generate covariance matrices from the resulting 
data-vectors. We also generate a Gaussian field covariance 
matrix to compare with the more realistic non-Gaussian 
analyses. 

Using these covariance matrices, we generate Fisher ma- 
trices for the £l m — as plane and show that the la, two- 
parameter contour for the non-linear Fisher matrix has an 
area that is 5.1 times larger than the theoretical prediction in 



both the 2-D and tomographic analyses. This indicates that 
the theoretical prediction is significantly under-predicting 
the errors on these parameters (the Figure-of-Merit is over- 
optimistic), even though their marginal errors are similar. 

To quantify if the contours generated using the Fisher 
matrix are a reasonable estimate of the true errors, we com- 
pared these outputs with a combined maximum likelihood 
analysis of the full simulation suite. The resulting contours 
in both the 2-D and tomographic analyses closely matched 
the contours generated with the simulation Fisher matrix. 
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From this, we can conclude that it is sensible to use a Fisher 
matrix analysis for parameter estimates, using a covariance 
matrix with all off-diagonal terms included. 

Based on the success of the non-linear Fisher errors in 
matching the maximum likelihood estimates in the Sl m — 
plane, we also perform both 2-D and 3-bin tomographic 
analyses to generate multi-parameter Fisher matrices of 
Q m , os, h, n s , wo and w a . We compare the non-linear 
Fisher matrix with the Gaussian Fisher matrix under the 
assumption that the simulation Fisher matrices are provid- 
ing accurate error estimates. With the survey parameters 
used in this paper, the 2-D analysis finds marginal errors 
on the parameters that are so large that they provide no 
constraining power. 

The tomographic analyses find reasonable marginal er- 
rors and demonstrate the value of the additional informa- 
tion obtained from the tomography. However, the marginal 
errors of the Gaussian Fisher matrix are larger than the non- 
linear Fisher marginal errors for every parameter which is 
a counterintuitive result. We show that the projected con- 
tours from the Gaussian Fisher matrix are also larger than 
the non-linear Fisher matrix, even though the volume of 
the non- linear Fisher is 3.7 times larger than the Gaussian 
Fisher. This warns us that the projected Fisher contours 
can be misleading over complex multi-variate spaces and 
that larger 2-D contours do not necessarily indicate a larger 
error volume. 
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